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The Einstein relation, relating the steady state fluctuation properties to the linear 
response to a perturbation, is considered for steady states of stochastic models with 
a finite state space . We show how an Einstein relation always holds if the steady 
state satisfies detailed balance. More generally, we consider nonequilibrium steady 
states where detailed balance does not hold and show how a generalisation of the 
Einstein relation may be derived in certain cases. In particular, for the asymmetric 
simple exclusion process and a driven diffusive dimer model, the external perturba- 
tion creates and annihilates particles thus breaking the particle conservation of the 
unperturbed model. 
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I. INTRODUCTION 

A guiding principle of equilibrium statistical mechanics is provided by the fluctuation- 
dissipation theorem, which relates the fluctuation properties of a system in equilibrium to the 
response of the system to an external perturbation. Its simplest form - the zero-frequency 
limit - is the Einstein relation which relates linear transport coefficients to spontaneous 
equilibrium fluctuations. For example the mobility /i and diffusion constant A of a Brownian 
particle interacting with its environment are related by 

// = /?A , (1) 

where (3 is the inverse temperature. This is a deep result with many implications. It enables 
one to reduce the calculation of response and transport coefficients (mobility, susceptibil- 
ity, etc.) to a zero- field calculation. Further, one can formally consider that spontaneous 
equilibrium fluctuations are caused by fictitious random forces. In particular the Einstein 
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relation dictates the forms of dissipative terms and noise in a phenomenological description 
of, for example, dynamic critical behaviour. 

A standard derivation of the Einstein relation (in the above form) comes from the 
Langevin equation describing Brownian motion. Here, the relation implies that diffusion is a 
direct consequence of the fluctuations of a Brownian particle, and, conversely, the existence 
of a frictional drag on the Brownian particle implies spontaneous equilibrium fluctuations of 
a random force. Because of the general applicability of the Einstein relation, one can make 
such statements for any transport problem for example, the existence of an electrical 
impedance implies spontaneous fluctuations of the voltage difference across its terminals P| . 

The most general formulations, applicable to quantum and classical models of equilibrium 
statistical mechanics, are due to Kubo 0, 4,1^, and Kadanoff and Martin ^,0]. The former 
relies on formal perturbation arguments for the linear response of a system near equilibrium; 
the latter approach is based on correlation functions. These approaches have been employed 
to derive the linearised hydrodynamics for the KLS model - a driven diffusive system and 
nonequilibrium generalisation of the Ising model - where the steady state satisfies detailed 
balance and the applied perturbation is the drive 0. 

Some attempts have been made to generalise the fluctuation-dissipation theorem to mod- 
els without detailed balance 0, 0, [l^ . These consider the relationship between dissipation 
and spontaneous steady state fluctuations and have found most success for steady states 
where the system can be divided into cells small on the macroscale, in each of which the 
dynamics satisfy detailed balance with respect to a local fixed density product measure. The 
global steady state can then be written in terms of these local product measure states. Such 
steady states are relevant when boundary driving breaks the global detailed balance of an 
equilibrium model. Also, a fluctuation theorem for nonequilibrium steady states has been 
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However, it is not known in general how to derive an Einstein relation for steady states 
of driven nonequilibrium systems where detailed balance is not satisfied. Since no such 
guiding principle exists in nonequilibrium statistical mechanics, we seek here to generalise the 
Einstein relation to nonequilibrium steady states of lattice-based stochastic models governed 
by a Master equation. 

In what follows, we review in section |H] the approach introduced by Derrida that 
yields general expressions for a current and diffusion constant for systems described by 
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a Master equation. We consider stochastic models with a finite state space and use this 
approach to investigate the existence of Einstein relations. In section UTTl we consider steady 
states where detailed balance is satisfied and show generally that an Einstein relation holds. 
In section HVl we consider nonequilibrium steady states where detailed balance is not obeyed 
and discuss the conditions under which an Einstein relation holds. We show that the usual 
Einstein relation holds for a simple class of models that includes boundary-driven symmetric 
exclusion (diffusing particles with hard-core interaction). Beyond this class, we consider 
bulk-driven models such as the asymmetric simple exclusion process (ASEP) and driven 
dimer diffusion. In particular, we obtain an Einstein relation for the ASEP that relates 
the diffusion constant to the response of the current to an applied field that creates and 
annihilates particles. We conclude in section IVl 



II. CALCULATION OF THE CURRENT AND DIFFUSION CONSTANT 



In this section, we outline how to obtain general expressions for a current J and diffusion 
constant A, for stochastic models defined on a finite state space. The approach is given in 



detail by Derrida et. al. 



m 



15|; here we just present the key equations. 



We consider a random variable Yt which increases and decreases according to specified 
dynamical events. For example, for a system of random walkers, if we take Yt as the total 
displacement, then Yt increases when a walker steps to the right and decreases when one steps 
to the left. In the long-time limit, the mean and variance of Yt increase linearly in time, 
and the constants of proportionality are interpreted respectively as a current (sometimes 
referred to as a velocity) and a diffusion constant. This is always true for finite systems, 
regardless of whether or not the infinite system exhibits anomalous diffusion i.e. when the 
mean or variance of Yt increases nonlinearly in time. 

The moments of Yt are obtained from Pt{C,Y), which is the probability that at time t 
the system is in a configuration C and the random variable takes the value Yt = Y. The 
time evolution of Pt{C, Y) is governed by a Master equation 

f^Pt{C,Y) = J2 [{MoiC,C')PtiC',Y) + M,iC,C')PtiC',Y-l) + M_i(C, C')Pt(C', F+1)} 
c 

- {MoiC',C)PtiC,Y) + M^iC',C)PtiC,Y) + M_,iC',C)PtiC,Y)}] , (2) 
where Mo(C, C) is a transition rate from the configuration C to the configuration C without 
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contributing to Yf, Mi{C,C') contains transitions from C to C which increase Yf by one 
and M_i{C,C') contains transitions from C to C which decrease Yt by one. The diagonal 
elements Mo{C,C), Mi{C,C) and M_i(C,C) are defined to be zero. It is straightforward to 
obtain evolution equations for the first two moments of Yt, (Yt) and (Y^), by multiplying 
both sides of (j2I) by y and Y'^. The angled brackets here denote an average over steady 
state initial conditions and all histories of the dynamics. (The formalism can be extended to 

n 

include other initial conditions il6|.) These evolution equations are expressed in a compact 
fashion in terms of PtiC) = YIiy ^tiC^ ^) ^"^^ ItiC) = J2y ^^tiC^ VtiC) is the probability 
of being in configuration C at time t; qt{C)/pt{C) is the average value of Yt given the system 
is in configuration C at time t. These quantities tend toward the asymptotic steady-state 
forms |3| 

Pt{C) ^ p{C) and qt{C) ^ Jtp{C) + r{C) , (3) 
in the long-time limit. 

Here J is a current defined by J = dt(Yt) and is found, by substituting the asymptotics 
into the evolution equation for (Yt), to be 

J = Y, Wi{C,C') - M_,{C,C)]p{C') . (4) 

C,C' 

The diffusion constant A associated with the current J is defined as A = dt{{Y^) — (Yt)"^). 
Then substituting the asymptotics Q into the evolution equations for (Yt) and (Y^) yields 

A = ^ [Ml (C, C) + M_i (C, C')] p{C') + 2 5^ [Ml (C, C) - M_i (C, C')] r (C) - 2 J r (C) . 

C,C' C,C' c 

(5) 

The quantities p{C) and r(C) are determined by substituting the asymptotic forms dHJ 
into the evolution equations for pt{C) and qt{C), which are obtained from 0: p{C) satisfies 

J2 [M{C, C')p{C') - M(C', C)p{C)] = , (6) 
c 

where the total transition matrix M = Mq + Mi + M_i, and r(C) satisfies 

5^[M(C,C')r(C')-M(C',C)r(C)] = Jp(C)-^[Mi(C,C')-M_i(C,C')]p(C') • (7) 
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Although (|7j) only fixes r{C) up to a term proportional to p{C), this term will cancel from 
dHl), so we take without loss of generality 



5^r(C)=0. 



(8) 



III. EINSTEIN RELATION FOR SYSTEMS SATISFYING DETAILED 

BALANCE 



We now illustrate how one may exploit the approach outlined in Section |n] in order to 
obtain an Einstein relation in the case where the dynamics obey detailed balance. 
Firstly, the condition under which detailed balance holds is 

M{C,C')p{C') = M{C',C)p{C). (9) 

Under the assumption that if Yf increases by one on going from C to C, then Yf decreases 
by one on going from C to C, (jH)) implies 

Mi(C, C')p(C') = M_i(C', C)p(C) , (10) 
M_i(C, C')p{C') = M,{C', C)p{C) . (11) 

Using these in equation Q implies that the current J = 0. Therefore, from (|7j), we find 
that r(C) satisfies 

J2 [MiC, C')t{C') - M(C', C)r(C)] = - J] [M,(C, C) - M_,{C, C')] p{C') . (12) 
c c 

Next , we consider a perturbation which couples only to dynamics which evolve the random 
variable Yt i.e. the perturbed transition matrix is written 



M(C, C) = Mo(C, C) + e^Mi(C, C) + e-^M_i(C, C) 



(13) 



where 7, which measures the strength of the perturbation, is small (and where M[C, C) = 0). 
Then, to first order in 7, the steady-state Master equation ^ requires that 



E 



A,(C,C')^-M(C',C)*"'' 



^7 



(97 



^7 



^7 



(14) 



We can use (tT3|l in the rhs of this equation along with the detailed balance conditions (jlOpilll 
to write 



E 



M(C,C')^-A/(C',C)*"'' 



(97 



^7 



-2 5^ [Ml (C, C) - M.i (C, C')] p(C') . (15) 



6 



By comparing equations (fT^ and (fT3jl we deduce 

2KC) ^ Mil . (16) 

Therefore we can express the diffusion constant, given by (0), in terms of the response of 
the steady state to the apphed field: 

A = J2 [Ml {C,C') + M_^{C,C')]p{C') + J2 [Ml {C, C) - M^^{C, C')] . (17) 

C,C' C,C' 

From the expression of the response dJjd^ of the current to the perturbation ()13|). to 
first order in 7, is the same as the rhs of (fT7|) . Thus, 

which is the usual equilibrium Einstein relation. 

Thus we have demonstrated an Einstein relation for all finite-state stochastic models 
which obey detailed balance (in the absence of the perturbing field). Previously Einstein 
relations have been derived for particular models in this class. For example, the symmetric 



exclusion process 



17| and the zero-field repton model of polymer dynamics|l 



IV. EINSTEIN RELATION FOR NONEQUILIBRIUM STEADY STATES 

In the absence of detailed balance, the structure of the steady state is not known in 
general and it is not known how to formulate an Einstein relation. 

For the case of detailed balance we showed in the previous section that the Einstein 
relation resulted from the simple relationship given in equation (fT^ . between r(C) and the 
response of p(C) to the applied field. Inspired by this, in the following in cases where detailed 
balance doesn't hold, our aim is to find a perturbation, parameterised by a rate 7, such that 

2KC) = . (19) 

where is a constant of proportionality. 
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A. Simple case of usual Einstein relation 

For a perturbation of the form (fT^ . it is clear by comparing ((Tj) and (fT^ that (fT^ holds 
for the class of models which satisfies the condition 

J = Y,[M^{C\C) - M^r{C' ,C)] , (20) 
c 

for all configurations C. Using the condition (jH}, the diffusion constant for this class of 
models can be written 

A = J] [Mi(C,C') + M_i(C,C')]p(C') , (21) 

which using the expression for J (jH) yields the usual Einstein relation (fTKj) . Thus pUjl is a 
simple condition for the usual Einstein relation to hold for nonequilibrium systems. 

We now show that the condition ()20|) holds for the boundary driven symmetric exclusion 
process. This is a stochastic model defined on a lattice of sites, where each site i is either 
occupied by a particle (indicated by the variable rii taking value 1) or vacant (indicated by 
nj = 0). In one dimension, the dynamics in the bulk are defined such that particles hop to 
the right (which we represent as the process 1 ^ 1) or the left (the process 01 10) 
with rate 1. At the left-hand boundary site, particles are injected with rate 1 — and 
removed with rate pi and at the right-hand boundary site, particles are removed with rate 
Pr and injected with rate 1 — pr. Thus the system is in contact with boundary reservoirs 
at densities pl at the left-hand boundary and pr at the right-hand boundary. This forces a 
current through the system (provided Pl 7^ Pr) and in the steady state the density profile 
is linear. 

We define Yt such that it is increased (decreased) every time a particle is added (removed) 
at the left-hand boundary, every time a particle hops to the right (left), or every time a 
particle is removed (added) at the right-hand boundary. In the absence of the perturbation, 
the exact current is [1^ 

J = PR- PL. (22) 
(Note that the current across a bond is J/(iV + 1).) To verify that piU) holds, this is to be 
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compared with 

[Mi(C', C) - M_,{C', C)] =(1 - - ^i) - PlUi + Pru^ - (1 - Pr){1 - n^) 

c 

+ Y K(l - 'f^i+i) - (1 - • (23) 

For any configuration C = {rii}, the terms in the sum over i cancel except near the bound- 
aries, and one obtains 

J2 [MiiC, C) - M.,{C', C)] =PR-PL, (24) 
c 

for every configuration C. Therefore the condition ^J]^ holds and the Einstein relation (fTH|) 
follows. 

The class of models for which ()20|1 holds includes the boundary driven symmetric zero- 
range process and other models where detailed balance is broken by the boundary dynamics 
i.e. models which would satisfy detailed balance if periodic boundary conditions were im- 
posed. 

B. Einstein relation for the ASEP 

We now consider the asymmetric simple exclusion process (ASEP). In one dimension the 
model, with periodic boundary conditions, is defined such that particles hop only to the 
right with rate 1 (i.e. the process 10^01). In the steady state, all configurations occur 
with equal likelihood, so p{C) is given by the binomial coefficient, 

-1 

= ^N,M ' (25) 

for a ring of sites containing M particles. For the ASEP neither detailed balance nor 
the condition ()2()|1 are satisfied. However, by considering a perturbation which creates and 
annihilates particles, an Einstein relation can be derived in the following way. 

We define Yt such that it is incremented every time any particle performs a hop, so the 
current (0} is given by 

(26) 
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and using equation (|7j), r(C) must satisfy 

J2 [MiC, C')r(C') - M(C', C)r (C)] = iNZ^\,Z^_,^M-i - ho) p{C) , (27) 



where /lo = ~ ^i+i) is the number of 1 configurations across all bonds in the 

system configuration C. 

Now consider the perturbed transition matrix for the ASEP, 

M(C,C') =Mi(C,C') + M^(C,C') , (28) 

where the elements of Mi(C,C') describe the unperturbed ASEP and the elements of 
My(C,C') describe the following processes that occur at all pairs of nearest neighbour sites 

11^10, 11^01, 

00^10, 00^01, (29) 

where 7 is a rate and {cj} measure the relative strengths of the four processes. Since this 
perturbation creates and annihilates particles, the different particle sectors (i.e. configura- 
tions with the same number of particles) are now connected by the dynamics. This modifies 
the steady state: if we define P{M) to be the (normalised) probability that the system is in 
the sector containing M particles, then, for 7 — > 0, the probability that the system is in the 
configuration C in the M-particle sector, pm{C), can be written as 

Pm{C) = P{M)p{C) (30) 

where p{C) is the steady state probability in the absence of the perturbation. P{M) is 
determined by the following balance condition with respect to transitions between particle 
sectors in the limit 7 — > 0: 

(ei + 63) (1 1)m+i P{M + 1) = (e2 + £4) (0 )mP(M) , (31) 

where {■)m represents a correlation function calculated with the respect to the steady state 
f)25|) in the M-particle sector. Since in the limit 7 — all configurations within each 
sector are equally likely, these correlation functions are easily found to be (ll)Af+i = 
^n\i+i^n-2,m-^ and (0 0)^./ = ZJj^j^Zn-2,m- Consequently the solution of (p?T|) is read- 
ily obtained as 

P{M) = A^i ^-^-^ Zn,m Zm-2,m-i , (32) 
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leading to the modified steady state puiC) given by 

Pm{C) = I ^ , ^ I AJj^Zn-2,m-i ) (33) 

where A^r is a normalisation, fixed by the requirement that YIm P{M) = 1. 
For the ASEP, (USD is shown to hold as follows. By (fHll . 



A/(C,C')%^-A/(C',C)*#> 

a7 a7 



(ei + e3)lupM{C) + (62 + e4)looPM{C) 



- eihoPM+i{C) - e^loiPM+iiC) 

- e2lioPM^i{C) - eJoiPM~i{C) . (34) 

The rhs is simplified by noting that Iiq = /qi and In + Iiq = M and /qo + ^lo = ^ ~ M, 
using (jnSD, and exploiting the identity Z^^m = Zn-2,m~'2 + 2^Ar--2,M-i + Zm-2,m- Thus, it 
is straightforward to show that the rhs of ()34|) is proportional to the rhs of equation ()27p. 
therefore r(C) is proportional to dpM{C)/d'j and the constant of proportionality Q is 

^ = ^nZ],' . (35) 



Hence the diffusion constant (0) can be expressed in terms of the response of the steady 
state to the perturbation, leading to an Einstein relation of the form 

8 J 

A = J + n— . (36) 

Thus for the ASEP we have shown that an Einstein relation (jHBP holds which expresses 
the diffusion constant as the current J plus the response of the current to perturbations (jSHD 
that create and annihilate particles. However it should be noted that not all perturbations 
that create and annihilate particles lead to For example creating and annihilating 

particles at a site regardless of the occupation of neighbouring sites does not satisfy (USD- 

It can be shown that a relation of the form ()19p holds for several variations of the ASEP 
under suitable perturbations which we now describe. 

a) Partially asymmetric exclusion: in the partially asymmetric exclusion process 
(PASEP) particles hop to the left with rate q and to the right with rate p (if the 
target sites are empty). This model can be treated, as above, by introducing the 
perturbation (jSHD) leading to an Einstein relation of the form 

A = ^J+{p~q)n^, (37) 
p — q Oj 

where the current J = {p — q)Z]^\jZ]\f_2^M-i and ^2 is given by (jSSD- 
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b) PASEP in higher dimension: an Einstein relation holds for the PASEP in any di- 

mension, on a hypercubic lattice with fully periodic boundary conditions. To illustrate 
this, consider a two dimensional square lattice on which particles hop up and to the 
right with rate p, and down and to the left with rate q. The random variable Yt is 
defined such that Yj increases whenever a particle hops up or to the right and decreases 
whenever a particle hops down or to the left. By considering a perturbation 

11-^10, 11^01, 00 10, 00^01, 

applied to horizontal pairs of nearest neighbour sites, and 

1^0 1_L^1 0^0 0^1 

11' 10' 1' 00' 

applied to vertical pairs of nearest neighbour sites, an Einstein relation of the form 
can be derived, where fl is given hj Q = K]^Z^^^. 

Further, if the dynamics are such that particles hop up and down with the same rate 
1, and to the right with rate p and to the left with rate g, then, again, the Einstein 
relation ()37|) can be derived. This is achieved by defining Yt such that Yt is increased 
when a particle hops to the right and decreased when a particle hops to the left (i.e. 
Yt is unaltered by hops up or down). The perturbation ()29|1 is then applied to pairs 
of nearest neighbour sites only along the axis of asymmetric hopping. 

c) A marked bond: if Yt is incremented only whenever a particle hops across a single, 

specified bond (the 'marked' bond), then a perturbation applied across this bond only 
and defined by the processes 

11^10, 00^10, 01-^10, (38) 

yields a relation of the form with Vt = AnZ]^'^j^. 

Finally we note that for sjTiimetric hopping, which satisfies detailed balance, the kinds of 
perturbations given in this section (i.e. those linking particle sectors without coupling to the 
original dynamics) do not satisfy the condition - rather, the appropriate perturbation 
couples only to the hopping dynamics, as discussed in section IIIII 
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C. Einstein relation for dimer diffusion 



In a similar way to the ASEP discussed in the previous subsection, we are also able to 
obtain an Einstein relation for (reconstituting) dimer diffusion: the process 

110 — ^011, (39) 

which occurs at all nearest neighbour triples of sites with rate 1. We assume the density 
is greater than one half so that there are no configurations where all particles are isolated. 
Again, in the steady state of this model, all configurations are equally likely; we define 
Yt such that it increases by one every time the process takes place (anywhere on the 
lattice), hence the current is given by J = NZ]^\jZn_.^^m-2- The Einstein relation (jHBj) is 
derived for a perturbation defined by the processes 

111^110, 100^110, 101^110. (40) 

As before, the steady state distribution is modified to accommodate a balance condition 
with respect to transitions between particle sectors c.f. ()31|) . Using the same reasoning as 
that applied to the ASEP, r(C) is shown to be proportional the the response of the steady 
state to the perturbation, where the constant of proportionality Q is 

Q = AnZ^^_^^j^j_^Z^^j^ , (41) 

and the diffusion constant and the response of the current to the perturbation are related 
by the Einstein relation given in equation ()36|) . 



D. Einstein relation for the boundary driven ASEP 

The boundary driven asymmetric simple exclusion process also satisfies an Einstein re- 
lation. In this model, particles can hop only to the right, and the multiple occupancy of 
any site is still forbidden. At the left-hand boundary site particles are injected with rate 
a and at the right-hand boundary site they are removed with rate j3. In this case, Yt is 
incremented every time a particle is added at the left-hand boundary site. Along the special 
line a + (3 = 1, all configurations in the same particle sector are equally likely in the steady 
state, and it is known that jl^ 

r{C) = a{l - a)^ , (42) 
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is a solution of ((Tj). This perturbation corresponds to a small increase in a away from the 
line a + /3 = 1, thereby changing slightly the density of the left-hand boundary reservoir. 
Hence, one finds an Einstein relation of the form 

Qj 

A = (2a- l)J + 2a(l -a)— . (43) 



V. CONCLUSION 



In this work we have investigated the existence of Einstein relations for stochastic, lattice 
models. Within the framework reviewed in section |n] we have shown that an Einstein 
relation will result if one can express the quantities r(C) as the response of the steady state 
probabilities p(C) to a perturbation. For the case of detailed balance we showed that this 
can always be done through the perturbation (fT!?|) which yields the usual Einstein relation 
fll8|) . Further for a class of nonequilibrium steady states where ()2U|) holds one again has the 
usual Einstein relation. This class includes models where detailed balance is broken only by 
..e bo.,*.. dy„..,. .c. . bou„d„y-dHven excision fiQ. 

Turning to systems where detailed balance is lacking even under periodic boundary condi- 
tions we have found nonequilibrium generalisations of the Einstein relation for some specific 
models. In these cases the perturbation creates and annihilates particles, breaking the con- 
servation law of the unperturbed dynamics. 

The models for which we succeeded in finding an Einstein relation had the simplifying 
property that in each particle sector all configurations are equally likely. It remains a 
challenge to establish whether Einstein relations hold for more general nonequilibrium steady 
states. Natural starting points would be steady states that have a simple structure such as 
the one dimensional KLS model (if. 
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